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We clarify the behavior of curvature perturbations in a nonlinear theory in case the infiaton tem- 
porarily stops during inflation. We focus on the evolution of curvature perturbation on superhorizon 
scales by adopting the spatial gradient expansion and show that the nonlinear theory, called the be- 
yond <5iV-formalism for a general single scalar field as the next-leading order in the expansion. Both 
the leading-order in the expansion (<5iV-formalism) and our nonlinear theory include the solutions 
of full-nonlinear orders in the standard perturbative expansion. Additionally, in our formalism, we 
can deal with the time evolution in contrast to <5iV-formalism, where curvature perturbations remain 
just constant, and show decaying modes do not couple with growing modes as similar to the case 
with linear theory. We can conclude that although the decaying mode diverges when (j> vanishes, 
there appears no trouble for both the linear and nonlinear theory since these modes will vanish at 
late times. 

PACS numbers: 98.80.-k, 98.90.Cq 



I. INTRODUCTION 

Recent observations of the cosmic microwave back- 
ground anisotropy [l| show very good agreement of the 
observational data with the prediction of standard infla- 
tionary cosmology, that is, adiabatic Gaussian random 
primordial fluctuations with an almost scale-invariant 
spectrum generated from quantum fluctuations of an in- 
fiaton field during inflation 2 8] . The amplitude of cur- 
vature perturbation on a comoving slicing, 1Z C , is given 
by the formula 7Z C « H 2 /\<fr\ evaluated at time of hori- 
zon crossing t = ifc when the wavenumber k satisfying 
k — aH, where <p and H are the infiaton and the Hubble 
parameter during inflation, respectively. The reason why 
it gives an almost scale-invariant spectrum is that both 
H and <p change very slowly during slow-roll infiaton. 

The main purpose of this paper is to clarify what hap- 
pens when the infiaton stops during that, namely, <f> « 0, 
using a nonlinear perturbation theory. Such a situation 
naturally occurs in oscillating inflation or the chaotic new 
inflation models [9l-[lTj]. For example, it has been shown 
by Damour and Mukhanov that oscillating inflation is 
realized as the infiaton oscillates around a minimum of 
a nonconvex potential 0] . In another example of chaotic 
new inflation fill ], it has been pointed that the infiaton 
changes its direction of motion if model parameters are 
approximately chosen. 

In both examples, if we use the above formula for 
the amplitude of primordial curvature perturbation, it 
apparently diverges when 6 vanishes. However, Seto, 
Yokoyama and Kodama [12J have shown that even in 
the case that slow-roll conditions arc violated, if a new 
formula is applied to this case, the amplitude has still 
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finite value, where the time derivative of the scalar field 
is replaced by the potential gradient given as 1Z C oc 
3H 3 /V'(4>) at t = tfc. In this study, they have inves- 
tigated the evolution of curvature perturbations in the 
linear theory and shown a decaying mode can diverge at 
and around temporary stopping of the infiaton fl2i Tfl3|. 
However, since the decaying mode vanishes at sufficient 
late times, there appears no trouble in the linear pertur- 
bation theory. 

When wc take nonlinear effects into account, the de- 
caying modes can couple with the growing modes in gen- 
eral to convert into growing modes through such effects. 
Therefore the ill-behavior of the decaying mode may leave 
an observable trace if nonlinear perturbation is incorpo- 
rated. The purpose of this paper is to clarify the behavior 
of curvature perturbations in a nonlinear theory in case 
the infiaton temporarily stops during inflation. 

In order to incorporate nonlincarity of curvature per- 
turbation, we focus on the evolution on superhorizon 
scales and consider a nonlinear cosmological perturbation 
theory by adopting a gradient expansion approach [3]. 
As for the leading-order in the expansion, <5iV-formalism 
is a powerful tool to calculate the nonlincarity 
of primordial curvature perturbatio ns ( recently much at- 
tention to as their non-Gaussianity [19() since it includes 
the solutions of full-nonlinear orders in the standard per- 
turbative expansion, but this is just a lowest-order and 
in this formalism, we should ignore all decaying modes. 
Therefore, we have to use the next-leading order in the 
expansion, which was recently formulated by one of us, 
the so-called beyond 5N- formalism [2(|. In our formal- 
ism, which we will briefly review in the following section, 
there exists decaying and growing modes, having their 
time-dependences and they lead to time variations of su- 
perhorizon curvature perturbations. We will show such 
decaying mode also diverges, but they will vanish due to 
inflationary expansion in the same way as in the linear 
theory, when (f> vanishes. 
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The rest of the paper is organized as follows. In Sec. 
II, we review the full-nonlinear cosmological perturbation 
theory of superhorizon curvature perturbations. Then 
we discuss the growing and decaying modes in both the 
linear and nonlinear theories in Sec. Ill and discuss what 
happens on temporary stopping of the inflaton in Sec. IV. 
Section V is devoted to the conclusion. 



II. BEYOND 57V-FORMALISM 

In this section, we will briefly review the nonlinear the- 
ory of cosmological perturbations valid up to 0(e 2 ) in the 
spatial gradient expansion and follow the previous works 
pol [2lj , where e is the ratio of the Hubble length scale 
1/H to the characteristic length scale of perturbations L, 
used as a small expansion parameter, e = 1/(HL), of the 
superhorizon scales. First of all, we show the main result 
in our formula for the nonlinear curvature perturbation, 
7? NL 

K" h " + 2^K" h ' + ^KW[Tl™ L ] = 0(e 4 ), (2.1) 

which shows two full-nonlinear effects; 

1. Nonlinear variable: including full- nonlinear 
curvature perturbation, SN 

2. Source term: K^[R^ h ] is a nonlinear function of 
curvature perturbations. 

In (|2.1[) , the prime denotes conformal time derivative and 
z is a well-known Mukhanov-Sasaki variable which will 
be seen later as (|2.24[) . The explicit forms of both the 
definition of 7Z^ h and the source term K^ 2 '[X], that is 
the Ricci scalar of the metric X, will be also seen later, in 
(|2.23|) and in (|2.30p , respectively. Of course, in the linear 
limit, it can be reduced to the well-known equation for 
the curvature perturbation on comoving hypersurfaces 

^Lin" + 2 -Tl^ n/ ~ c 2 s A[R% in ] = . (2.2) 

z 

We will briefly summarize our formula and show the 
above results in the following. Throughout this paper we 
consider a minimally-coupled single scalar field described 
by an action of the form 

1= J d i x^^P(X,(j)), (2.3) 

where X = — g^ d^(j)d u 4>. Note that we do not assume 
the explicit forms of both kinetic term and its potential, 
that can be given as arbitrary function of P(X, </>). 

We adopt the ADM decomposition and employ the gra- 
dient expansion. In the ADM decomposition, the metric 
is expressed as 



where a is the lapse function, ft is the shift vector and 
Latin indices run over 1,2,3. The equations of motion 
corresponding to a and ft lead to constraint equations. 
Components of the spatial metric 7^ are dynamical vari- 
ables and the corresponding equations of motion are re- 
duced to a set of first-order differential equations with 
respect to the time t. We introduce the extrinsic curva- 
ture Kij defined by 

Kij =~ ($7« - Aft - Aft) , (2-5) 

where D is the covariant derivative compatible with the 
spatial metric 7^ . As a result, the basic equations are re- 
duced to the first-order equations for the dynamical vari- 
ables (jijjKij), with the two constraint equations (the 
so-called Hamiltonian and Momentum constraint). We 
further decompose the spatial metric and the extrinsic 
curvature as 

2 ic ~ 

7ij = a " e 7*3) 

Kij = a 2 e 2 < (^Kjij + I y J , (2.6) 

where a(t) is the scale factor of the background FRW 
universe and det^- = 1. 

Next, we will employ the gradient expansion. In this 
approach we introduce a flat FRW universe (a(t), 4>o{t)) 
as a background. As discussed in the first part of this 
section, we consider the perturbations on superhorizon 
scales, that is, L is longer than the Hubble length scale 
1/H of the background, i.e. HL ^ 1. Therefore, we 
consider e = 1/(HL) as a small expansion parameter and 
systematically expand our equations by e, considering a 
spatial derivative acted on perturbations is of order O(e). 

We assume the condition for the gradient expansion: 

9t% = 0(e 2 ). (2.7) 

This corresponds to assuming the absence of any decay- 
ing modes at the leading-order in the expansion, namely, 
the absence of spatially homogeneous anisotropy. This 
is justified in most of the inflationary models in which 
the number of e-folds of inflation N is much larger than 
the number required to solve the horizon and flatness 
problem, N ^> 60. This assumption is sufficient to allow 
us discuss behavior of decaying modes when the inflaton 
stops, since all time dependent solutions at the leading 
order are reduced to just decaying modes and there exists 
no observable trace at late times. 

When we focus on a contribution arising from the 
scalar-type perturbations, we may choose the gauge in 
which 7jj approaches the flat metric, 

% (t 00) = 5ij, (2.8) 

where in reality the limit t — > 00 may be reasonably 
interpreted as an epoch close to the end of inflation. We 
take the comoving slicing, time- orthogonal gauge: 



ds 2 = -a 2 dt 2 + ^{dx 1 + p i dt)(dx j + P j dt), (2.4) 



5«MM*)=ft(i,**)=0(e 3 ), (2.9) 
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where 5cj> = 4> — 4>o denotes a fluctuation of a scalar field. 
The subscript c denotes this gauge throughout this paper. 

Now we turn to the problem of properly defining a non- 
linear curvature perturbation to 0(e 2 ) accuracy. Here- 
after we will use the expression 1Z C on comoving slices to 
denote it. Let us consider the linear curvature perturba- 
tion which is given as 



n 



Lin 



rrLin 

11 rp 



Y. 



(2.10) 



where, following the notation in |23j . the spatial metric 
in the linear limit is expressed as 

7v = « 2 (*u + 2#L in ^<% + 2#T ir %-) , (2-11) 

with Y being scalar harmonics with eigenvalue k 2 in 
Fourier space satisfying 



(A + k 2 )Y = 0. 



and 



Yij — k 



didj - -<%A 



Y 



(2.12) 



(2.13) 



These expressions in the linear theory correspond to the 
metric components in our notation as 



C = H^Y, % 



(2.14) 



Notice that the variable £ c reduces to IZ];" 1 at leading- 
order in the gradient expansion, but not at second-order 
as f|2 . 10[) and it will be also similar to the nonlinear the- 
ory. 

Thus to define a nonlinear generalization of the linear 
curvature perturbation (|2.10p . we need nonlinear gener- 
alizations of H^Y and HtY . Our nonlinear £ is an ap- 
parent natural generalization of H^ m Y, 



h l y = c 



(2.15) 



As for HtY , however, the generalization is non-trivial. It 
corresponds to the 0(e 2 ) part of 7^ and we have obtained 
a general solution of the dynamical equation for 7^ as a 
first-order differential equation in [2(| [2l| and the time- 
dependent part includes the following solution; 



7ii (i) 3 C. 



(2) 



dt' 



a 3 (i')' 



(2.16) 



with the Momentum constraint: 



by using (12.17[) . As shown in [20(, it can be done by in- 
troducing the inverse Laplacian operator A -1 on the flat 
background and we defined the nonlinear generalization 
of HtY as 



H T Y = E = --A" 



<-)' 



- 3f( Ve 3 * <0, (ln7), 



.(2.18) 



It is easy to see that E 3 which we expected. 

At leading-order, the only non-trivial quantities for the 
spatial metric, £ and 7^ , are given by 



( = e (0 Hx k ) + O(e 2 ), 



and 



4V) + 0(e 2 ), 



(2.19) 



(2.20) 



where l^{x k ) is an arbitrary function of the spatial co- 
ordinates {x k } (k = 1,2,3) and f£\x k ) is a (3 x 3)- 
matrix function of the spatial coordinates with a unit 
determinant, respectively. Throughout this paper, this 
leading-order of spatial metric can be chosen as 



(o) 



(2.21) 



consistent with the gauge condition of (|2.8[) . On the 
other hand, £ (0) represents a conserved comoving cur- 
vature perturbation, equivalent to a fluctuation of the 
number of e-folds, which is denoted by the so-called SN 
term from some final uniform density (or comoving) hy- 
persurface to the initial flat hypersurface at t = t*, 



e<® =6N(t„x l ) 



(2.22) 



With these definitions of H^Y and HtY, we can define 
the nonlinear curvature perturbation valid up through 
0(e 2 ) as 



n 



NL 



Cc + 



E r 



(2.23) 



It is easy to show that this nonlinear quantity can be re- 
duced to (|2.10[) in the linear limit. As clear from (|2.18l) . 
finding HtY generally requires a spatially non-local oper- 
ation, however, in the comoving slicing, time-orthogonal 
gauge with the asymptotic condition on the spatial coor- 
dinates (|2.8p . we find it is possible to obtain the explicit 
form of HtY without any non-local operation as seen in 

Next, we can derive a nonlinear second-order differ- 
ential equation that 1Z^ L (|2.23j) satisfies at 0(e 2 ) ac- 
curacy by introducing the conformal time 77, defined by 
drj = dt/a(t) and the Mukhanov-Sasaki variable [22], 



Q f(0) d J 



(2.17) 



The explicit forms of solutions can be seen in [21] . Here 
we attach the superscript (m) to a quantity of 0(e m ), 
and both and will be denoted as the leading- 
order metric in (|2.19p and (|2.20p . Our aim is to derive the 
scalar-type solution from the tensor cjj in (|2.16p 



a_ f p + P 
H 



(2.24) 



where notice that c s is the speed of sound for the gauge 
invariant scalar perturbation in the linear theory [24j, 
given by 



P 



x 



Px + 2P XX X 



(2.25) 
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where the subscript X represents derivative with respect 
to X. The result can be reduced to a simple equation of 
the form (|2.1|) as a natural extension of the linear ver- 
sion (|2.2p . We also obtain the solution of the nonlinear 
equation (|2.ip as 



eV> + l[F(rj)-F.]KW 

+ [D(rf) - D,]C<V +0(e 4 ), (2.26) 



where 



D(rj) = m 
F{rf) 



v * 2 (V 
drf r-n 



z 2 M , , 
—) d71 ' 

" z 2 c 2 s (v"W 



(2.27) 



<i-< 



Here D* = D(t]^,), F* = F(r]*) and denotes the con- 
formal Hubble parameter H = dlna/drj at r/ = r/» which 
we take the time as some after the horizon crossing. Note 
that t — > oo corresponds to i] — > in the conformal time. 
Thus the functions D and F vanish asymptotically at late 
times, D(0) = F(0) = 0. Deviation of the solution (|2T2"6"j) 
can be easily understood as follows. The second-order 
differential equation (|2.1[) contains two solutions (even 
though its independent relation appears only for the lin- 
ear theory), i.e. decaying mode and growing mode. We 
can find that the function D (rj) satisfies 



(2.28) 



D" + 2-D' = 0. 

z 



in the long- wavelength limit, i.e. no source term in (|2.1j) . 
It will be seen that it corresponds to the decaying mode 
in the linear theory in the next section. On the other 
hand, the function F(rf) corresponds to the source term 
in (|2.ip . satisfying 



F" + 2-F' ■ 



0. 



(2.29) 



as the 0(e 2 ) correction to a constant mode at the leading- 
order, i.e. as the growing mode in the linear theory, which 
is taken the form 1 + F(r/)A' (2 ) + (9(e 4 ). 

Moreover the equation (|2.1[) includes two 'constants' of 
integration, or arbitrary spatial functions, which in gen- 
eral appear as the initial conditions, namely, the initial 
value and its time derivative. Let us consider the spa- 
tial functions, which we have introduced as C^ 1 and 
A'( 2 '. Here the last one is related to the Ricci scalar of 
the Oth-order spatial metric as 



(2)rKo)i _ 



K (2) [e ( 



R 



e Oi 



= -2(2A£<°) + <r-W (0) <9/ 0) )e^ (0) , (2.30) 



where we have used = Sij from (|2.21[) . Then we have 
the two arbitrary spatial functions: £( ' and C^ 2 \ which 
are related to the number of physical degrees of freedom 



for the initial conditions. Therefore £^ and corre- 
spond to the initial conditions determined by matching a 
solution of n-th order perturbation solved inside the hori- 
zon to this superhorizon solution at r/ = 77* . Notice that 
represents SN term as seen in (|2.22[) and orig- 
inally comes from the decaying mode of the fluctuation 
of the scalar field 1201 . 



III. GROWING AND DECAYING MODES 

In this section, firstly, let us consider the growing and 
decaying modes in the linear theory. In the linear theory, 
the curvature perturbation on comoving hypersurfaccs 
follows (|2.2[) . As usual, we consider it in Fourier space, 



K 



Lin' 



2—K L c 

z 



ik 2 n^ n = 0. 



(3.1) 



Real space expressions (|2.2p will be recovered by the re- 
placement k 2 — > — A. This equation has two independent 
solutions, conventionally called a growing mode and a de- 
caying mode. 

The growing mode is a constant at the leading-order 
in the long- wavelength approximation or equivalently the 
spatial gradient expansion. Then in terms of the growing 
mode solution u, the decaying mode solution v can be 
given as [25[ 



, z 2 Mu 2 M 
^ z 2 {r ] ')u 2 {7 1 ') [6 ' 



Note that this expression is correct for any order in the 
gradient expansion in the linear theory. 

The general solution of a curvature perturbation is 
written in terms of their linear combinations as 



(3.3) 



where the coefficients a Lm and (3 Lm may be assumed to 
satisfy a Lm + /3 Lln = 1 without loss of generality. Note 
that the assumption of the gradient expansion (|2.7[) cor- 
responds to the condition, 



1 



0(e 2 ). 



(3.4) 



This means, as mentioned before, that the decaying mode 
at leading-order in the gradient expansion has already 
decayed after horizon crossing. 

Therefore the decaying mode solutions can be auto- 
matically obtained as following (|3.2I) . if we obtain the 
growing mode solutions. Let us solve for the growing 
mode solution. In accordance with the gradient expan- 
sion, we set 



(3.5) 



At the leading-order in the gradient expansion, the grow- 
ing mode solution u^ ' is just a constant. Then inserting 
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the above expansion with =const. to the equation 
of motion fl3.1[) gives iteratively 



-l +2— <+i 



~C 2 s U n 



(3.6) 



As shown in [25j], 0(k 2 ) corrections to can be written 



= u (°) 



cf > + C y z>D( V ) + k 2 F{rf) 



,(2) 



(3.7) 



where the integrals D(rj) and F(ji) have been given 
in (f2~27)) , satisfying ((235)) and (|2T2l)|) , respectively, as 

(2) (2) 

similar to the nonlinear theory, and C[ and C\ are 
constants of integration. We fix the two constants as 
C[ = C 2 = so that tr 2 ' is proportional to the inte- 
gral F(jj) at 0(k 2 ) accuracy 1 . Hence we find 



1^(77) = k 2 u ( ^F( V ) 



(3.8) 



As for the decaying mode, because of (|3.4[) we only 
need the leading-order solution. Since we may replace D 
with D in (|3.2[) . we immediately find 



,(o) 



(3.9) 



Thus from (|3.8p and (|3.9p . the general linear solution 
valid up through 0(e 2 ) is obtained as linear combination 
of constant u^°\ growing mode and decaying mode 
!)", which are proportional to F(rf) and D(rj), respec- 
tively 

As for the nonlinear theory of cosmological perturba- 
tions, the solution of (|2.26[) is also shown as growing and 
decaying modes, respectively. We can find that the func- 
tion D(r)) and F(ij) satisfy and (|2~2^| . respectively 
and they take the same forms, respectively as in the linear 
theory. Therefore we can interpret that they correspond 
to the decaying mode in the long- wavelength limit, and 
the growing mode taken the form 1 + k 2 F(t]) + 0(fc 4 ), 
where F(rj) is the k 2 correction to the growing (i.e., con- 
stant) mode, respectively. In our nonlinear theory, note 
that time derivative takes the same form as shown in (|2.1[) 
and (|2.2[) . hence the decaying mode can not couple with 
the growing mode as similar to the linear theory because 
of the method of gradient expansion (i.e. time deriva- 
tive takes as a linear operator). The difference from the 
linear theory is the source term , i.e. the Ricci scalar of 

the leading order metric K^ 2 \ which can be reduced to 
fc 2^Lin in Fourier 

space as the source term in the linear 

theory (|3~T|) . 



1 If wc fix the two arbitrary constants as cj 2 ' = and Cj = 
—k 2 F,/D, so that «(?7* ) = holds at 0(k 2 ) accuracy, it is the 
case of [25l in which they discussed an enhancement of curvature 
perturbation on superhorizon scales due to suddenly change of 
the infiaton potential's slope, and its nonlinear effect also can be 
studied by matching the linear solution of [25| to our nonlinear 
solution in 12011 



IV. CROSSING OF <j> = 

We consider the case when <p (or z) crosses zero in 
our nonlinear theory. The gradient expansion allows us 
to discuss in a similar way as the linear theory [l2|, [25j . 
For simplicity, we assume that z changes the sign only 
once at r\ = r] . Hereafter, we consider a canonical single 
scalar field, however, the same discussion can be done 
in the case of a non-canonical single scalar field, when 
P X X » 0. 

In the vicinity of r\ = rj , z can be expressed as 

z = 4iV-Vo), (4.1) 
where z' Q = 2/(770). Hence the equation for 1Z^ L becomes 
2 d 



d?_ 

drf rj — 770 dr\ 
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-K^['R^] + 0{e i ). 



(4.2) 



The two independent solutions in the linear theory (not 
guaranteed in the nonlinear theory) can be found as 



*(0) (Vi^ 2 >[7^ L ](77-77o) 2 + 



(4.3) 



„ „ C ^ (-L- l -K^[U^ ](„ -vo) + ■■■)■ (4.4) 

We consider u and v should be chosen as the growing 
mode and decaying mode, respectively, and u remains 
constant across the epoch 77 = 770. The second term in 
(|4.3[) can be obtained by the integral ^(77), which in this 
case is given by 



F(rj) oc lim (77 - r) y 



(4.5) 



and shown to be still well defined in the crossing of cj> = 0. 
The final value of the growing (or non-decaying) mode at 
late times will take a constant (i.e. SN term). 

The singularity will appear in the first term in (14.41) . 
arising from the integral D(rj). It can be expressed as in 
the case of linear theory and for 77 > 770, we obtain as 



f° dn' 1 

Jr, z 4 W" 



(4.6) 



We conclude that this term diverges in the limit 77 — s- 
770 + 0, however, this is just a decaying mode, then it will 
vanish definitely at late times as 



D{rf) (x a 3 -» 0, with 77 -> 0. 



(4.7) 



Hence we can see that no problem will occur for both the 
linear and nonlinear theory. 
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V. CONCLUDING REMARKS 

We clarify what happens when the inflaton stops dur- 
ing inflation for nonlinear cosmological perturbation the- 
ory. We focus on the evolution on the supcrhorizon scales 
and review our nonlinear theory, called the beyond 5N- 
formalism for a general single scalar field as the next- 
leading order in the gradient expansion. In our nonlinear 
theory, we can deal with the time evolution in contrast 
to (5./V-formalism where curvature perturbations remain 
just constant. 

As a summary of our formula, note that time derivative 
takes the same form as shown in (|2.1[) and (|2.2I) . hence 
the decaying mode can not couple with the growing mode 
as similar to the linear theory because of the method of 
gradient expansion, i.e. time derivative takes as a linear 



operator. The difference from the linear theory is the 
source term , i.e. the Ricci scalar of the leading-order 
metric K {2 \ which can be reduced to k 2 7i^ ln in Fourier 
space as the source term in the linear theory (|3.ip . 

We can conclude that although the decaying mode di- 
verges in the limit of time when <f) vanishes, there appears 
no trouble for both the linear and nonlinear theory since 
this mode will vanish definitely at late times. 



VI. ACKNOWLEDGEMENTS 

This work was supported by JSPS Grant-in-Aid for 
Young Scientists (B) No. 21740192. and for Scientific 
Research No. 19340054. 



[1] E. Komatsu et al. . larXiv: 1000 538 [astro-ph.CO]. 

[2] A. H. Guth, Phys. Rev. D 23 347 (1981). 

[3] A. A. Starobinsky, Phys. Lett. B 91 99 (1980). 

[4] K. Sato, Mon. Not. Roy. Astron. Soc. 195 467 (1981). 

[5] V. F. Mukhanov and G. V. Chibisov, JETP Lett. 33 532 

(1981) [Pisma Zh. Eksp. Teor. Fiz. 33 549 (1981)]. 
[6] S. W. Hawking, Phys. Lett. B 115 295 (1982). 
[7] A. A. Starobinsky, Phys. Lett. B 117 175 (1982). 
[8] A. H. Guth and S. Y. Pi, Phys. Rev. Lett. 49 1110 (1982). 
[9] T. Damour and V. F. Mukhanov, Phys. Rev. Lett. 80, 

3440 (1998). 

[10] A. R. Liddle and A. Mazumdar, Phys. Rev. D 58, 083508 
(1998); A. Taruya, Phys. Rev. D 59, 103505 (1999). 

[11] J. Yokoyama, Phys. Rev. D 58, 083510 (1998); 59, 
107303 (1999). 

[12] O. Seto, J. Yokoyama and H. Kodama, Phys. Rev. D 61, 
103504 (2000). 

[13] H. Kodama and T. Hamazaki, Prog. Theor. Phys. 96, 
949 (1996). 

[14] D. S. Salopek and J. R. Bond, Phys. Rev. D 42, 3936 
(1990). 

[15] A. A. Starobinsky, JETP Lett. 42, 152 (1985) [Pisma Zh. 



Eksp. Teor. Fiz. 42, 124 (1985)]. 
[16] Y. Nambu and A. Taruya, Class. Quant. Grav. 13, 705 
(1996). 

[17] M. Sasaki and E. D. Stewart, Prog. Theor. Phys. 95, 71 
(1996). 

[18] D. H. Lyth and Y. Rodriguez, Phys. Rev. Lett. 95, 
121302 (2005). 

[19] N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, 

Phys. Rept. 402, 103 (2004). 
[20] Y. Takamizu, S. Mukohyama, M. Sasaki and Y. Tanaka, 

JCAP 1006, 019 (2010). 
[21] Y. Takamizu and S. Mukohyama, JCAP 0901, 013 

(2009). 

[22] V. F. Mukhanov, H. A. Feldman and R. H. Branden- 

berger, Phys. Rept. 215, 203 (1992). 
[23] H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78 

1 (1984). 

[24] J. Garriga and V. F. Mukhanov, Phys. Lett. B 458, 219 
(1999). 

[25] S. M. Leach, M. Sasaki, D. Wands and A. R. Liddle, 
Phys. Rev. D 64, 023512 (2001). 



